Decoherence vs entanglement in coined quantum walks 
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Quantum versions of random walks on the line and cycle show a quadratic improvement 
in their spreading rate and mixing times respectively. The addition of decoherence to the 
" quantum walk produces a more uniform distribution on the line, and even faster mixing on 

the cycle by removing the need for time-averaging to obtain a uniform distribution. We 
calculate numerically the entanglement between the coin and the position of the quantum 
walker and show that the optimal decoherence rates are such that all the entanglement is 



(N 



just removed by the time the final measurement is made. 
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I. INTRODUCTION AND BACKGROUND 

Simple quantum generalisations of classical random walks spread quadratically faster on the 
line QjSj) an d m i x quadratically faster on the cycle jjj. These promising examples of a quantum 
speed up were soon followed by several algorithms based on quantum walks. Shenvi et al. 0] 
proved a quantum walk can solve the unsorted database search problem quadratically faster, and 
Childs et al. proved an exponential speed up for crossing a particular type of graph. Several 
more algorithms have followed these, Ambainis |6j gives an overview of quantum walk algorithms, 
and Kempe Q] provides an introductory review of quantum walks and their properties. Physical 
implementations of quantum walks have also been proposed , as sensitive tests of coherent 

control over quantum particles. Optical implementation dates back to Bouwmeester et al. fill ] . 



which can be viewed as a quantum walk on the line with many photons 12|, [13|]. A quantum walk 
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on a cycle with four nodes has been implemented on a NMR (nuclear magnetic resonance) quantum 
computer using three qubits, two for the position label (in binary) and one for the coin 14j |. 

Quantum walks on simple one-dimensional structures remain a fertile testing ground for further 
research. It was shown numerically by Kendon and Tregenna [15] that the addition of decoherence 
or measurements to the quantum walk dynamics can optimise the spreading and mixing properties 
for quantum walks on both the line and the cycle (recently proved by Richter la, ll7|). A detailed 
survey of the effects of decoherence in quantum walks can be found in la]. In this work we 
look closely at how the interplay between quantum evolution and decoherence or measurements 
produces optimal computational properties. In particular, we calculate the entanglement between 
the coin and the position of the quantum walker and observe how it varies as the quantum walk 
evolution unfolds. This study is restricted to quantum walks taking place in discrete time and 
space, using a quantum coin to control the choice of direction. Thus we can use the entanglement 
between the coin and position as an indication of how "quantum" the walk is as it progresses. 
A different method for assessing "quantumness" would be needed for continuous-time quantum 
walks, as introduced for quantum algorithms by Farhi and Gutmann [l9| . 

In the following subsections we describe the background material: simple quantum walks, the 
decoherence model we use, and the entanglement measure (negativity) we employed. The results 
of our study are then described in fTH first for the line, then for cycles, and finally a summary of 
the common elements of our findings. 



A. Simple, coined, quantum walks 

A discrete-time (coined) quantum walk dynamics consists of a quantum "coin toss" operation 
C, followed by a shift operation S to move the quantum walker to a new position. These are 
repeated alternately for T steps of the quantum walk, and the final position of the quantum walker 
is measured. For quantum walks on the line and the cycle, we have just two choices of which way 
to step, so the quantum coin is a two state system. We write \x, c) for a quantum walker at position 
x with a coin in state c G {+1, — !}• For a walk on the line, x £ 7* and for a walk on a cycle of 
size N, we have x £ Zjy. The full state of the quantum walk |^(t)) at time t can be written as a 
superposition of terms in each basis state \x, c), 

|*(t))=X>*.c(t)l*>C>, (1) 

x,c 

where ip x ,c(t) E C and the normalisation is c \ip XjC (t)\ 2 = 1. When the quantum walk is measured 
(in the basis just defined), the probability of finding the quantum walker at position y with the 
coin in state b is given by 

P(y,b,t) = \(y,b\*(t))\ 2 = \4>y,b(t)\ 2 - (2) 
The coin toss C and shift S are defined in terms of their action on the basis states \x,c), 

S\x, c) = \x + c, c) (3) 



C\x,c) = (\x, — c) + c\x, c))/v2 



(4) 
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One can add more general bias or phase into the coin toss operation, see, for example, 20], but 
this does not greatly change the basic properties of the quantum walk on a line or cycle, so we will 
consider only the unbiased case in this paper. For a quantum walk starting at position x = with 
the coin in a superposition state (|— 1) + i|+l))/y / 2, (where i = \J — 1) we can write a quantum 
walk of T steps as 



*(T)) = (SC) J {|0, -1) + i\0, +1)} /V2 (5) 



The solution for |^(T)) may be obtained by various methods such as Fourier analysis [2| and path 
counting [l| , and has been studied extensively. The key result is that spreading on the line proceeds 
linearly with the number of time steps. We can use the standard deviation <tq(T) of the probability 
distribution to quantify the spreading rate. For a quantum walk on the line, 



a Q (T) = ^2x 2 P(x,c,T)^i [ l-—j T, (6) 



asymptotically in the limit of large T [2]. In contrast, for the classical random walk, the standard 
deviation is (Jc(T) = \[T . 

On the cycle, we are interested in mixing times rather than spreading. Mixing times can be 
defined in a number of different ways, we choose the main definition given in yj], 

M(e) = min {T | V t > T : \\P(x, t) - P u \\ tv < e} (7) 

where P u is the limiting distribution over the cycle, and the total variational distance (TVD) is 
defined as 

\\P(x,T) - P u \\ tv = \P(x,T) - P u \. (8) 

X 

A classical random walk on the cycle mixes to within e of the uniform distribution in time propor- 
tional to iV 2 log(l/e), where e can be chosen arbitrarily small. 

Pure quantum walks, on the other hand, do not mix to a stationary distribution. Their deter- 
ministic dynamics ensures they continue to oscillate indefinitely. There are several ways to obtain 
mixing behaviour, first explored by Aharonov et al. [3|. By defining a time-averaged probability 
distribution for the quantum walk, 

T-l 

P(x,c,T) = ^P(x,c,i) (9) 

t=0 



they proved that P(x, c, T) does converge to a stationary distribution on a cycle, and that on odd- 
sized cycles the stationary distribution is uniform. A mixing time can be defined for P(x,c,T), 

M(e) = min{r | V t > T : \\P(x,t) - P u \\ tv < e} . (10) 

and Aharonov et al. 0] proved that, for odd-sized cycles, M(e) is bounded above by 0(e~ 3 N log N), 



almost quadratically faster (in N) than a classical random walk. Kendon and Tregenna 151 ] . 
observed numerically that M(e) ~ O(Nfe), this has been recently confirmed analytically by Richter 
0. 
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FIG. 1: Numerical data for a quantum walk on the line of 100 steps, with decoherence applied to both the 
coin and position, for various decoherence rates between zero (quantum) and one (classical) as shown in the 
key. 



Notice that we pay a price for the time-averaging: the scaling with the precision e is now linear 
instead of logarithmic. Aharonov et al. Q provide a fix for this in the form of a "warm start". 
The quantum walk is run several times, each repetition starting from the final state of the previous 
run. A small number of such repetitions is sufficient to reduce the scaling of the mixing time M(e) 
to logarithmic in e. 

Both the quantum walk on the line and the cycle thus provide a quadratic speed up over 
classical random walks. This quadratic speed up does not carry over to all quantum walks on 
higher dimensional structures, see, for example, 
ubiquitous this behaviour really is. 



21 



22, 



23l . |24|. It remains an open question how 



B. Decoherence in quantum walks 

We will consider decoherence in the form of randomly-occurring uncorrelated non-unitary events 
added to the quantum walk dynamics already described. The evolution of the quantum walk must 
now be described using a density operator p(t) given by 

p(i + l) = (l-pJSCpCtjC+St+pJ^PySCp^C+Stpt. (11) 

j 

Here Fj is a projection that represents the action of the non-unitary decoherence events and p is the 
probability of applying the decoherence per time step, or, completely equivalent mathematically, 
to a weak coupling between the quantum walk system and a Markovian environment with coupling 
strength p. For a pure state, the density operator p(t) = \^{t)){^{t)\, it thus has the normalisation 
Tr[p(t)] = 1. For p = 0, equation (jlljl reduces to the pure quantum walk described in §1 A\ and 
for p = 1, to a classical random walk. 

Full analytical solution of a non-unitary quan tum walk has been done only for a few special cases. 



For quantum walks on the line, Brun et al. 25f] analysed the case of random measurements on the 



coin only. While analytical solution is challenging, equation (jlip lends itself readily to numerical 
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FIG. 2: Numerical data for the time-averaged mixing time M e (p) on cycles of size N = 48, N — 49, N = 50 
and N = 51, for coin (dotted), position (dashed) and both (solid) subject to decoherence, using e = 0.002. 



simulation since p, S and C can be manipulated as complex matrices, and the Pj generally remove 
some or all of the off-diagonal entries in p. Kendon and Tregenna 1_5J] evolved equation (jllH 
numerically for various choices of ¥j\ projection onto the position space, projection onto the coin 
space, and projection of both coin and position. In all cases, the spreading rate is reduced, in the 
long time limit 25], it becomes proportional to vf instead of proportional to T. More interesting 
behaviour is seen for intermediate times and decoherence rates p with decoherence applied to the 
position, or to both the position and coin. Kendon and Tregenna observed that, for 2 < pT < 5, 
the distribution becomes very close to uniform while retaining the full quantum linear spreading 
rate, see figured! With decoherence applied to the coin only, the distribution retains a cusp shape 



151 ]. To quantify the distance from a uniform distribution, the TVD given by equation ([8]) can 



be used, this time with P u defined to be a top-hat of appropriate width x £ {±T/\/2}, see[l|, Q]. 
Since the quantum walk (in common with the classical random walk on which it is based) has 
the property that at odd(even) time steps the walker will only be found at odd(even)-numbered 
locations, the top-hat definition of P u also incorporates this property. 

For a quantum walk on the cycle subjected to decoherence, the mixing behaviour is dramatically 
improved provided the decoherence is applied to the position. Decoherence guarantees mixing to 
the uniform distribution, and a similar judicious choice of 2 < pN < 5 produces a minimum 



mixing time 151 ] . well below the classical value. Decoherence applied only to the coin does cause 



the quantum walk on a cycle to mix, but not significantly faster than a classical random walk. 
Furthermore, time-averaging is no longer necessary, mixing occurs in time 0(N log(l/e). This has 
recently been proved by Richter |1.7]. This mixing behaviour is illustrated in figures [2] and [3] for 
cycles of size 48, 49, 50 and 51. Figure [2] shows time-averaged mixing times while figure [3] shows 
mixing times without any time averaging. Note the scales are different between the two figures, 
the mixing times are longer when time-averaging is applied, reflecting the change from logarithmic 
to linear dependence on 1/e. 
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FIG. 3: Numerical data for the mixing time, without time-averaging, M e (p) on cycles of size N = 48 
to N = 51, for decoherence applied to coin only (dotted), position only (dashed) and both (solid), using 
e = 0.002. 

C. Entanglement in mixed states 

Quantum walks with decoherence on the line and the cycle thus provide two examples where the 
optimal computational properties are obtained for a judicious combination of quantum dynamics 
and decoherence. In order to investigate this phenomenon in more detail, we asked how "quantum" 
the system is by the end of the quantum walk when the optimal amount of decoherence is applied. 
To quantify the "quantumness" , we used the entanglement between the coin and the position, 
which depends on the quantum correlations in the system. 



We chose our entanglement measure to be the negativity [26|, |27|, [28|] because this can be 
calculated numerically in a fairly straightforward manner for density operators such as p(t), and 
there are few options that meet this criterion. First we must choose a division of our system 
into two (or more) subsystems between which to identify the entanglement. For our quantum 
walk, the natural division is between the coin and the quantum walker's position. We note that 
the entanglement across this division will be the same whether we regard the quantum walk as a 
physical system with a qubit coin and a unary position, or as an algorithm running on a quantum 
computer with a single qubit for the coin and the position encoded in binary in the remainder of 
the quantum register. We perform a partial transpose on one subsystem to obtain a new matrix 
p'(t). For example, the partial transpose with respect to the coin subsystem is 

Pxcybit) = Pxb,yc(t) (12) 

where x, y are position indices and c, b are coin state indices. Next, we determine the spectrum of 
p'(t), denoted by {A^}. The normalisation of p{t) is carried over to p'(t), so ^ X' i = 1, but unlike 



p{t), it is possible for p'(t) to have negative eigenvalues. The negativity is defined [26|, |27|, [28| as 

E = \{^2\K\-^\, (13) 
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FIG. 4: Numerical data showing the negativity between the coin and position for decoherence rates p = 
(red), p = 0.05 (blue) and p = 0.1 (green) for a quantum walk on the line with decoherence on both coin 
and position, run for 40 time steps. 



which is just the sum of the negative eigenvalues. The negativity ranges between zero and one, 
with any non-zero value indicating entanglement is present. If the negativity is zero, it means the 
state is probably not entangled, but there can be exceptions 291 ] . The exceptions are known to be 
relatively rare in the set of all possible states 30J, and the entanglement they contain is difficult 
to apply to useful quantum tasks 29J]. For this study, we will not need the fine-grained detail of 
these possible exceptions. 

The entanglement in a pure state quantum walk has been studied previously, see for example, 



3l|. It fluctuates with each step, and eventually settles down to an asymptotic value that depends 



on the initial state of the quantum coin, and on any bias in the quantum coin operator C. The 
addition of decoherence smooths out this behaviour, see figure HI and steadily reduces the level of 
entanglement between the coin and the position. 



II. RESULTS AND CONCLUSIONS 



Our simulations simply take equation (jlip and evolve it numerically, calculating the negativity 
and, using an appropriate top hat or uniform distribution, the TVD for various types of decoherence 
and values of the decoherence rate p. We studied quantum walks on both the line and the cycle, 
for various initial states, lengths of walk, and sizes of cycles, with decoherence applied to the coin 
only, the position only, and to both the coin and position. These simulations can be accomplished 
on a desktop computer using straightforward code written in a language such as C++. We used 
various Apple Mac G4 and G5 computers and the versions of the GNU C compiler that come with 
the operating system OSX 10.4 and associated Developer Tools. Our results are summarised in 
the following subsections. 
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decoherence rate (per unit time step) 



FIG. 5: Numerical data showing the negativity (dotted, open symbols) and TVD (solid, filled symbols) 
for decoherence applied to the position only (blue triangles), the coin only (red circles) and to both coin 
and position (green squares) for a quantum walk on the line of 100 steps. Some negativity data is missing 
(missing symbols). 



A. Entanglement in the walk on the line 



Guided by the results in [15J] and [3l[ , we focused on the entanglement at the end of the quantum 
walk on the line, just before a final measurement would be made to find out where the quantum 
walker has ended up. We considered how the entanglement varies as the decoherence rate p is 
varied. As well as testing three cases of decoherence, applying it to both the coin and position, 
and also separately to just the coin or the position, we also tested two different starting states for 



the coin, both of which produce symmetrical distributions [22l. I32I]. these being (| — 1) + i|+l))/y2 
and cos(7r/8)|— 1) + sin(7r/8)|+l). The two initial states achieve a symmetrical distribution in 
fundamentally different ways, the former by combining two asymmetric distributions that do not 
interact, and the latter by arranging the interference exactly the right way to produce a symmetrical 
outcome. The distributions are not identical, and we observed that the TVD and negativity differ 
slightly, but not in any significant way. The negativity calculations require diagonalisation of 
large matrices (to find the eigenvalues), and we found for larger matrices (more time steps) our 



numerical routines (based on methods from Numerical Recipes [33() did not always manage to do 
this. Failure was always reported by the routines, and we have simply left gaps in the figures where 
the negativity data is missing. We obtained sufficient results to show the overall trends, as can be 
seen in a typical example from our results shown in figure Both the TVD from the optimal top 
hat distribution, and the negativity are plotted. The minimum in the TVD indicates the optimal 
decoherence rate p. For a quantum walk of 100 steps the optimal p is in the range 0.025 < p < 0.05. 
Note that the minimum for measurements on the coin only is much shallower. The distribution in 
this case retains a cusp shape [15j . The negativity is also higher than when decoherence is applied 
to the position of the walker, indicating that coin decoherence is less effective at removing the 
quantum correlations. 

If decoherence is applied to the position, with or without decoherence on the coin as well, the 
negativity drops to zero at p ~ 0.055, shortly after the optimal decoherence rate is reached. The 
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time step 



FIG. 6: Numerical data showing the negativity (dashed) and TVD (solid) for decoherence applied to the 
position for a quantum walk on the cycle of size 29 with near-optimal decoherence rates p = 0.2239 (blue 
triangles), p = 0.2511 (black circles) and p — 0.2818 (red squares). The pink horizontal line is at 1/29. 

optimal amount of decoherence is thus just about the amount required to remove the the quantum 
correlations from the system. 

B. Entanglement in the walk on a cycle 

For the mixing time on cycles, there are extra considerations. The mixing time depends not 
only on the size of the cycle, but also on how close to uniform one sets the threshold e, see equations 
(|7|) and (fTUl As already noted, pure quantum walks don't mix unless something is done to disrupt 
the pure quantum evolution. Both random decoherence and regular repeated measurements can 
efficiently change the behaviour into that of fast mixing to the uniform distribution [171 ] . Again 
decoherence applied to the coin only does not provide a significant improvement in the behaviour 
compared with decoherence applied to the position. Note also that the optimal rate of measurement 
2 < pN < 5 is independent of the threshold e, even though the main effect is to provide logarithmic 
scaling of the mixing time with e. 

We studied how the entanglement varies during the quantum walk on a cycle, with the deco- 
herence rate p chosen to be near-optimal. As with the walk on the line, we examined the three 
cases of decoherence applied to both the position and coin, and applied to just the position or coin 
separately. The results for a typical example, a cycle of size N = 29 with decoherence applied to 
the position only, are shown in figured! While the actual mixing time is determined by the choice 
of e, we have indicated the position of 1/N in figure El this being the probability of finding the 
walker at one location in a uniform distribution. The time at which the TVD drops below this line 
is around the time the entanglement also drops to zero. The variability of the TVD shows that the 
mixing time is not a smooth function of e, so we cannot expect to determine a more precise result. 
Nonetheless, we have the same qualitative behaviour as we found for the quantum walk on a line: 
the optimal amount of decoherence is just enough to remove all of the quantum correlations. 
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C. Discussion and conclusion 

For the purpose of uniform sampling, the optimal quantum walk on the line and the cycle is a 
carefully balanced combination of quantum dynamics with decoherence or measurements providing 
randomness during the evolution of the walk. We have shown that, while the quantum correlations 
are necessary to obtain linear spreading and mixing times, they must be neutralised to produce a 
uniform final distribution. This make intuitive sense: quantum correlations distort the distribution, 
giving it peaks and troughs, especially at the ends of the top hat [l|, Q]. On the other hand, if 
the decoherence rate is turned up until the classical random walk is obtained, classical correlations 
build up to produce the binomial distribution in which the quantum walker is more likely to be 
found nearer the starting point of the walk. A uniform distribution, by definition, has neither 
quantum nor classical correlations. 

The fact that decoherence applied only to the coin does not produce a good top hat distribution 
shows that this result is non-trivial. There is no guarantee that decoherence can be arranged to 
remove the quantum correlations almost completely before any significant classical correlations 
build up, but for the quantum walk on the line and the cycle it is possible to do this by applying 
decoherence to the position. When decoherence is applied to the coin only, classical correlations 
have started to build up before all the quantum correlations are removed. 

Applied to quantum computation, this result has intriguing implications. For the task of uniform 
sampling, the optimal way to use a quantum walk is not a pure quantum process, it is a slightly 
noisy process which has lost all entanglement by the end of the computation. Given the difficulty 
of arranging for perfect noiseless quantum operations, it might be easier to implement an optimal 
quantum walk since it doesn't have to be perfect. It also begs the question whether there are 
other computational tasks for which a noisy quantum process is optimal. In our search for perfect 
quantum computation we may have overlooked the possibility that some useful tasks may not need 
it. A communications task in which added randomness improves the performance is already known: 
distilling a secret key from imperfect quantum exchanges 0, 35]. Here the added noise reduces 



the information the eavesdropper can obtain about the secret key, and leads to better key rates 
than a deterministic key distillation process. 
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